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Abstract 

The robustness properties of bipartite entanglement in systems of iV bosons distributed 
in M different modes are analyzed using a definition of separability based on commut- 
ing algebras of observables, a natural choice when dealing with identical particles. 
Within this framework, expressions for the robustness and generalized robustness of 
entanglement can be explicitly given for large classes of boson states: their entangle- 
ment content results in general much more stable than that of distinguishable particles 
states. Using these results, the geometrical structure of the space of N boson states 
can be explicitly addressed. 



1 Introduction 

When dealing with many-body systems made of identical particles, the usual definitions of 
separability and entanglement appear problematic since the natural particle tensor product 
structure on which these notions are based is no longer available. This comes from the fact 
that in such systems the microscopic constituents can not be singly addressed, nor their 
properties directly measured [HE]- 

This observation points to the need of generalized notions of separability and entangle- 
ment not explicitly referring to the set of system states, or more in general to the "particle" 
aspect of first quantization: they should rather be based on the second quantized description 
of many-body systems, in terms of the algebra of observables and to the behavior of the as- 
sociated correlation functions. This "dual" point of view stems from the fact that in systems 
of identical particles there is not a preferred notion of separability, it can be meaningful only 
in reference to a given choice of (commuting) sets of observables [3]-[E]- 
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This new approach to separability and entanglement have been formalized in [9]- [12]: it 
is valid in all situations, while reducing to the standard one for systems of distinguishable 
particles^] in particular, it has been successfully applied to the treatment of the behavior of 
trapped ultracold bosonic gases, giving rise to new, testable prediction in quantum metrology 

0-HB. 

As in the case of systems of distinguishable particles]! suitable criteria able to detect 
non-classical correlations through the implementations of practical tests are needed in order 
to easily identify entangled bosonic many-body states [28j-[30j. In the case of bipartite 
entanglement, the operation of partial transposition [31]- [21] turns out to be again very 
useful; actually, it has been found that in general this operation gives rise to a much more 
exhaustive criterion for detecting bipartite entanglement than in the case of distinguishable 
particles [HI [12]. As a byproduct, this allows a rather complete classification of the structure 
of bipartite entangled states in systems composed of N bosons that can occupy M different 
modes [12], becoming completely exhaustive in some relevant special cases, as for M — 2. 

In the following we shall further explore the properties of bipartite entanglement in bosonic 
systems made of a fixed number of elementary constituents. We shall first study to what 
extent an entangled bosonic state results robust against its mixing with another state (sepa- 
rable or not): we shall find that in general, bosonic entanglement is much more robust than 
the one of distinguishable particles. In particular, we shall give an explicit expression for 
the so-called "robustness" [35J and upper bounds for the "generalized robustness" [36J. As 
a byproduct, a characterization of the geometry of the space of bosonic states will also be 
given; the structure of this space results much richer than in the case of systems with iV 
distinguishable constituents. One of the most striking results is that the totally mixed sepa- 
rable state, proportional to the unit matrix, does no longer lay in the interior of the subspace 
of separable states: in any of its neighborhood entangled states can always be found. 



2 Entanglement of multimode boson systems 

As mentioned above, we shall focus on bosonic many-body systems made of ./V elementary 
constituents that can occupy M different modes. From the physical point of view, this 
is a quite general framework, relevant in the study of quite different systems in quantum 
optics, atom and condensed matter physics. For instance, this theoretical paradigm is of 
special importance in modelling the behavior of ultracold bosonic gases confined in multi- 
site optical lattices, that are becoming so relevant in the study of quantum phase transitions 
and other quantum many-body effects (e.g., see [251 [2B], [3TJ-[44J and references therein). 
In order to properly describe the N boson system, let us thus introduce creation a\ and 

lr The notion of entanglement in many-body systems has been widely discussed in the recent literature 
{e.g., see P2]-[2I]); nevertheless, we stress that only a limited part of those results are relevant for identical 
particles systems. 

2 For general reviews on the role of quantum correlations in systems with large number of constituents 
see H3-[22]- 
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annihilation operators dj, % = 1,2,..., M, for the M different modes that the bosons can 
occupy, obeying the standard canonical commutation relations, [a*, a]] = 6{j. The total 
Hilbert space H. of the system is then spanned by the many-body Fock states, obtained by 
applying creation operators to the vacuum: 

|ni,n 2 , . . . ,n M ) = ) (at)" 1 (4)" 2 ■ ■ ■ (4)"" |0> , (1) 

the integers rii,n 2 , . . . ,n M representing the occupation numbers of the different modes. Since 
the number of bosons is fixed, the total number operator a l a « i s a conserved quantity 
and the occupation numbers must satisfy the additional constraint Yli=i n i = in other 
words, all states must contain exactly N particles. As a consequence, the dimension D of 
the system Hilbert space H is finite; one easily finds: D = { ■ 

In addition, the set of polynomials in all creation and annihilation operators, {a\, a^}, 
i = 1,2, . . . , M, form an algebra that, together with its norm-closure, coincides with the 
algebra B(H) of bounded operators]! the observables of the systems are part of this algebra. 

When dealing with systems of identical particles, instead of focusing on partitions of 
the Hilbert space H, it seems natural to define the notion of bipartite entanglement by 
the presence of non-classical correlations among observables belonging to two commuting 
subalgebras Ai and Ai of B(H) [9]. Quite in general, one can then introduce the following 
definition: 

Definition 1. An algebraic bipartition of the algebra 13(1-1) is any pair (A\,A2) of 
commuting subalgebras of 13(H), A\,Ai C 13(H). 

More explicitly, a bipartition of the M-oscillator algebra 13(H) can be given by splitting the 
collection of creation and annihilation operators into two disjoint sets, {af , Oj \i — 1, 2 . . . , m} 
and {a], aj, \ j — m + 1, m + 2, . . . , M}; it is thus uniquely determined by the choice of the 
integer m, with < m < M. All polynomials in the first set (together with their norm- 
closures) form a subalgebra Ai, while the remaining set analogously generates a subalgebra 
Ai- Since operators pertaining to different modes commute, one sees that any element of 
the subalgebra A\ commutes with each element of A%, in short [Ai.,.4.2] = 0. 

Remark 1: i) Note that there is no loss of generality in assuming the modes forming the 
subalgebras A± (and A2) to be contiguous; if in the chosen bipartition this is not the case, 
one can always re-label the modes in such a way to achieve this convenient ordering. 
ii) Further, when the two commuting algebras Ai and A2 are generated only by a subset 
M' < M of modes, one can simply proceed as if the iV boson system would contain just 
the used M' modes, since all operators in B(H) pertaining to the modes not involved in 
the bipartition commute with any element of the two subalgebras A\ and A2, and therefore 
effectively act as "spectators" . As a consequence, all the considerations and results discussed 

3 The algebra B(H) is generated by the so-called Weyl operators; all polynomials in the in the creation 
and annihilation operators are obtained from them by proper differentiation |45l 146] . 
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below holds also in this situation, provided one replaces the total number of modes M with 
M', the actual number of modes used in the chosen bipartition. □ 

Having introduced the notion of algebraic bipartition (Ai,A 2 ) of the system operator 
algebra £>("H), one can now define the notion of local observable: 

Definition 2. An element (operator) of B(l-C) is said to be (Ai, A2) -local, i.e. local with 
respect to a given bipartition (A%, A2), if it is the product A 1 A 2 of an element A x of A\ and 
another A2 in A 2 . 

From this notion of operator locality, a natural definition of state separability and entangle- 
ment follows [9]: 

Definition 3. A state p (density matrix) will be called separable with respect to the 
bipartition (Ai,A 2 ), in short (Ai, A2) -separable, if the expectation Tr[pv4iA 2 ] of any local 
operator A\A 2 can be decomposed into a linear convex combination of products of expecta- 
tions: 

Tr [p A 1 A 2 ] = h Tr [p® Aj Tr [p k 2) A 2 ] , \ k > , ^ \ k = 1 , (2) 

k k 

where } and {pjP } are collections of admissible states for the systems; otherwise the 
state p is said to be entangled with respect the bipartition (A\, A 2 ). 

Remark 2: i) This generalized definition of separability can easily be extended to the case 
of more than two partitions; for instance, in the case of an n-partition, Eq.(J5]) above would 
extend to: 

Tr[pA 1 A 2 ---A n ]=^A fc Tr[pi 1) A 1 ]Tr[pl 2) A 2 ]---Tr^ n) A n ], X k > , ^ X k = 1 . (3) 

k k 

ii) When dealing with systems of distinguishable particles, one finds that Definition 3 gives 
the standard notion of separability 0. 

in) In this respect, it should be noticed that when dealing with systems of identical particles, 
there is no a priori given, natural partition, so that questions about entanglement and 
separability, non-locality and locality are meaningful only with reference to a specific choice 
of commuting algebraic sets of observables [3]- [12]; this general observation, often overlooked 
in the literature, is at the basis of the definitions given in (T2J) and ([3]). 

iv) A special situations is represented by pure states [12]. In fact, when dealing with pure 
states, instead of general statistical mixtures, and bipartitions that involve the whole algebra 
£>("H), the separability condition in ([2]) (and similarly for ([3])) simplify, becoming: 

Tr[pA 1 A 2 ] = Tr[pW A 1 ] Tr[p^ A 2 ] , p = |^)<VI , (4) 

with p^\ p^ projections on the restrictions of \ip) to the first, respectively second, partition; 
in other terms, separable, pure states turn out to be just product states. □ 
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Examples of N bosons pure separable states are the Fock states. Using the notation and 
specifications introduced before (cf. ([!])), they can be recast in the form: 



m M 

\ki,...,k m ;k m+1 ,...,k M ), ^ki = k, 2J kj — N — k , <k < N , (5) 

i=l j=m+l 

where k indicates the number of bosons in the first partition; by varying it together with 
the integers ki, these states generate the whole Hilbert space "H. This basis states can be 
relabeled in a different, more convenient way as: 

fk + m-l\ , fN — k + M — m — l\ 

\k,a;N-k,a') , a = l,2,..., , , o' = 1, 2, . . . , AT , ,(6) 



k J' ---••^ N _ k 

where, as before, the integer k represents the number of particles found in the first m modes, 
while a counts the different ways in which those particles can fill those modes; similarly, 
a' labels the ways in which the remaining N — k particles can occupy the other M — in 
modesQ In this new labelling, the property of orthonormality of the states in ([5]) becomes: 
(k, a;N - k, a'\l, t;N -I, r') = 5 M 5 aT 8 a < T >. 

For fixed k, the basis vectors {\k, a; N — k, a')} span a subspace "H^ of dimension D k D N _ k , 
where for later convenience we have defined (cf. (JH]) above): 

Remark 3: Note that the space 1-Lk is naturally isomorphic to the tensor product space 
C Dfe ® C DiV - fe ; through this isomorphism, the states \k,a;N — k,a') can then be identified 
with the corresponding basis states of the form \k, a) £g> | iV — k, a'). □ 

By summing over all values of k, thanks to a known binomial summation formula [47], one 
recovers the dimension D of the whole Hilbert space "H : 

Y,D k D N _ k = D=i ) . (8) 

fc=o ^ ' 

Using this notation, a generic mixed state p can then be written as: 

N N 

P=^2 Pkaa',irr' \k,a;N - k,a')(l,r;N -l,r'\ , ^ ^ Pw.w = 1 • (9) 

k,l=0 a;a',T,T' k=0 <j,a' 

In general, to determine whether a given density matrix p can be written in separable form 
is a hard task and one is forced to rely on suitable separability tests. One of the most useful 



4 Clearly, we need two extra labels a and a' for each value of A;, so that these labels (as well the range of 
values they take) are in general fc-dependent: in order to keep the notation as a simple as possible, in the 
following these dependences will be tacitly understood. 
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such criteria involves the operation of partial transposition [311 E2]: a state p for which the 
partially transposed density matrix p is no longer positive is surely entangled. This lack of 
positivity can be witnessed by the so-called negativity [3B1 SSI EE] : 

M(P)= - Tr[p]) , HplU = Tr [v^5] • (10) 
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which is nonvanishing only in presence of a non positive p. Although this criterion is not 
exhaustive (there are entangled states that remain positive under partial transposition), 
it results much more reliable in systems made of identical particles [TT| IT2] , Indeed, the 
operation of partial transposition gives a necessary and sufficient criteria for entanglement 
detection for very general classes of bosonic states e.g. in presence of only two modes 
(M = 2), or, in the generic case of arbitrary M, when the (A\, ^4.2)-bipartition is such that 
the algebra A± is generated by creation and annihilation operators of just one mode, while 
the remaining M — 1 modes generates A2- 

Even more interestingly, it turns out that entangled iV-body bosonic states need to be of 
a definite, specific form [12] : 

Proposition 1. A generic (m, M — m)-mode bipartite state (0) is entangled if and only 
if it can not be cast in the following block diagonal form 

N N 

P = "^Pk Pk , J^Pfc = 1 , Tr[p fc ] = 1 , (11) 

k=0 k=0 

with 

Pk= Pkvv'Mr' \k,a;N- k,a')(k,r; N - k,r'\ , ^/w,w = 1 , (12) 



£7,(7' ,T,T" 



(i.e. at least one of its non-diagonal coefficients Pkaa' ,Itt' i_h 7^ I, J 's nonvanishing), or if it 
can, at least one of its diagonal blocks pk is non-separableu 

Proof. Assume first that the state p can not be written in block diagonal form; one can then 
show [12] that it is not left positive by the operation of partial transposition and therefore it 
is entangled. Next, take p in block diagonal form as in flTTJ), fTl2l) above. If all its blocks pk 
are separable, then clearly p itself results separable. Then, assume that at least one of the 
diagonal blocks is entangled. By mixing it with the remaining blocks as in ffTTl) will not spoil 
its entanglement since all blocks pk have support on orthogonal spaces; as a consequence, 
the state p results itself non-separable. □ 

Having found the general form of non-separable iV-boson states, one can next ask how 
robust is their entanglement content against mixture with other states. This question has 
been extensively studied for states of distinguishable particles ESI SHI EDI EH] ; in the next 
section we shall analyze to what extent the results obtained in that case can be extended to 
systems with a fixed number of bosons. 



5 For each block pk, separability is understood with reference to the isomorphic structure C Dk <E) C D " 
mentioned before (see Remark 3). 
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3 Robustness of entanglement 



Several measures of entanglement have been introduced in the literature with the aim of char- 
acterizing quantum correlations and its usefulness in specific applications [28]- [30]. Starting 
with the entanglement of formation, most of these measures point to the quantification of the 
entanglement content of a given state. A different approach to this general problem has been 
proposed in [351 EE] : the idea is to obtain information about the magnitude of non-classical 
correlations contained in a state p by studying how much it can be mixed with other states 
before becoming separable. 

More precisely, let us indicate with M. the set of all systems states and by S C M. that of 
separable ones; then, with reference to an arbitrary (m, M — m)-mode bipartition, one can 
introduce the following definition: 

Definition 4. Given a state p, its robustness of entanglement is defined by 

R(p) = infjt | t > 0, 3a eM' CM, for which 77 = jy^y e 5} , (13) 

i.e. it is the smallest, non-negative t such that a state a exists for which the (unnormalized) 
combination p + ta is separable. 

Actually, various forms of robustness have been introduced: they all share the definition 
(fl3|) . but differ in the choice of the subset A4' from which the mixing state a should be 
drawn. In particular, one talks of generalized robustness R g when a can be any state [36J, 
while simply of robustness R s when a must be separable [55] . 

All robustness defined in this way satisfy nice properties; more specifically, they result 
entanglement monotones, i.e. they are invariant under local operations and classical com- 
munication, and convex, R(\p\ + (1 — A)/?2j < A-R(pi) + (1 — A) R{p2)] further, R(p) in 
ffT3"j) is vanishing if and only if p itself is separable. Although implicitly proven in the case of 
states of distinguishable particles [33 [3B] , these properties hold true also in the case of N- 
boson systems. Nevertheless, there are striking differences in the behavior of the robustness 
of states of identical particles with respect to what is known in systems with distinguishable 
constituents. 

Let us first focus on the robustness R s (p), that measures how strong is the entanglement 
content of a state p when mixed with separable states. One finds: 

Proposition 2. The robustness of entanglement of a generic (m, M — m)-mode bipartite 
state p is given by 

N 

R s (p) = J2pk r M > ( 14 ) 

for states that are in block diagonal form as in OTj) . Q2j) . while it is infinitely large otherwise. 

Proof. From the results of Proposition 1, we know that separable iV-boson states must be 
block diagonal. If the state p is not in this form, it can never be made block diagonal by 
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mixing it with any separable one; therefore, in this case, the combination p + ta will never 
be separable, unless t is infinitely large. 

Next, consider the case in which the state p is in block diagonal form, i.e. it can be written 
as in ffTTT) . f[T2"l) . First, if p is separable, then clearly R s (p) = 0. For an entangled p, one 
can discuss each block p k separately: this is allowed since they have support on orthogonal 
Hilbert subspaces. Then, let us indicate by t k the robustness of block density matrix p k ; by 
Remark 3 and the definition of robustness, the numbers t k s are finite and positive, vanishing 
only when the corresponding state pk is separable [35]. More specifically, for each k, there 
exist separable states a k and rj k , such that: 

Pk + t k a k = (1 + 4)% • (15) 

Multiplying both sides of this relation by the positive number p k and then summing over k, 
one gets 

N 

p + tcr = (1 +t)r] , t = y^Pfctfc , (16) 

fc=0 

where the separable states a and r\ are explicitly given by 

N / + \ N /i i + 

E Pk^k \ / 1 + tfe 

fc=0 V 7 k=0 



Ok , r? = > J ( ) PkVk ■ (17) 



To prove that indeed t given in (TIB"]) is really the robustness of p, one needs to check that no 
better decomposition 

p + t'a' = (1 +tV , (18) 

with £' < t, exists. In order to show this, let proceed ad absurdum and assume that such 
decomposition can indeed be found. Since the states a' and rj' are separable, they must be 
block diagonal, i.e. of the form a' = J2 k qkCr' k and r( = Y J k r k r l'k wit h Y^k^k = J2 k r k = 1 
and a' k and r( k separable density matrices. By the orthogonality of the Hilbert subspaces 
with fixed k, from ( Fl8|) one then gets 

PkPk + t' q k a' k = (1 +t')r k r]' k , (19) 

and further, by taking its trace, p k + t' q k = (l+f) r k . In addition, from the previous identity, 
one sees that the combination 

Pk + t'^a' k , (20) 
Pk 

is separable. By definition of robustness of the block p k as given in f JT5|) . it then follows that: 

a - > h , (21) 

Pk 

ore equivalently, t k p k < t' q k . By summing over k, one then finds: t < t'\ this result 
is compatible with the initial assumption t > tf only if if coincides with t. Therefore, the 
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robustness of the block diagonal state p is indeed given by the weighted sum of the robustness 
of each block. □ 



The problem of finding the robustness R s {p) of a generic iV-boson state p is then reduced 
to the more manageable task of identifying the robustness of its diagonal blocks, which are 
finite-dimensional density matrices for which standard techniques and results can be used. 

Remark 4: i) A remarkable property of the robustness of entanglement of states describing 
distinguishable particles is that it is equal to the negativity for pure states. In the case of 
identical particles, this property does not hold anymore as the robustness of entanglement 
of non-block diagonal pure states is infinitely large. Nevertheless, one can easily show that 
in general for pure states: M{p) < R s (p)- 

ii) The robustness of entanglement of states that, with respect to the given partition, result 
mixtures of pure block states, p = ^2 k PkPk, Pk 2 = Pk, is equal to their negativity R s (p) = 
^2 k PkN(Pk) = A/"(p), since now R s (pk) = N(pk), as in the standard case (see [I2])- D 

More difficult appears the task of computing the generalized robustness R g (p) of a generic 
iV-boson state: only upper bounds can in general be given. In any case, note that in general 
Rg{p) — Rs{p) since the optimization procedure of Definition 4 is performed over a larger 
subset of states in the case of the generalized robustness. 

By fixing as before a (m, M — m)-mode bipartition, a first bound on R g (p) can be easily 
obtained. Let us extract from p its diagonal part po, as defined in terms of a Fock basis 
determined by the given bipartition (c/. fl§D), and call p^D = P — Pd the rest. By definition 
of separability, pr> is surely separable with respect to the chosen bipartition, so that an easy 
way to get a separable state by mixing p with another state a is to subtract from it its non- 
diagonal part. However, — p^n alone is not in general a density matrix since it might have 
negative eigenvalues. Let us denote by A the modulus of its largest negative eigenvalue; then 
the quantity Al — Pnd, where 1 is the identity matrix, will surely be positive and therefore, 
once normalized, can play the role of the density matrix a in the separable combination 
p + ta = pz> + \t. By taking the trace, one finds for the normalization factor t the following 
expression: t = XD, where D is as before the dimension of the total Hilbert space. By 
definition of robustness, it then follows that 

R g {p)<\D- (22) 

as a consequence, the generalized robustness of a generic iV-boson state is always finite. 
A different bound on R g (p) can be obtained using a refined decomposition for p, 

N N 

P = Pb+ Pnb , Pb = ^Pk Pk , ^Pk = 1 , Tr[p fc ] = 1 , (23) 

k=0 k=0 

where pB is the block diagonal part of p, whose blocks p k can be written as in f[T2"j) . while 
Pnb = P — Pb is the rest, containing the non block diagonal pieces. One can first ask for the 
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generalized robustness of p#, which is a bona fide state, being a normalized, positive matrix]^ 
Quite in general, one has: 

Proposition 3. The generalized robustness of entanglement of a generic (m, M — m)-mode 
bipartite state p given in block diagonal form as in ( TTT)) . ( fX2)) is given by 

N 

r 9 { P )=J2p* r m ■ ( 24 ) 

k=0 



Proof. It is the same as in Proposition 2; the only difference is that now the states are 
in general entangled. □ 

As a consequence, for a generic state as in (|23i) above, due to the presence of the additional 
term p NB , one surely has: R g (p) > R 9 (pb) = Y!k=oPk R g {pk)- 

To get an upper bound for R g (p), let us consider the form that the (unnormalized) density 
matrices a must take in order to make the combination p + o separable; by Definition 4, 
the generalized robustness of entanglement coincides with minimum value of their traces: 
R g (p) = inf { Tr [a] } . Since the combination p + o is separable, it must be in block diag- 
onal form; therefore, a must surely contain the contribution — Pnb- Further, it must also 
take care of the entanglement in the remaining block diagonal term ps] in view of the re- 
sults of Proposition 3, this can be obtained (and in an optimal way) by the contribution 
J2k=oPk Rg(.Pk) °~k, where cr^ is the optimal density matrix that makes the diagonal block pt 
separable. 

However, while J2k=oP k Rg(.Pk) a k is a positive matrix, in general p^B is not; therefore 
a should contain a further contribution, a (unnormalized) positive and separable matrix a 
curing the negativity induced by pnb- As a consequence, the generic form of the positive 
matrix a making the combination p + a separable is given by: 

N 

O- = ^2pkRg(pk)cTk- PNB + O- , (25) 
fc=0 

and the computation of the generalized robustness of p is reduced to the determination of 
the optimal 5"; indeed: 

N 

R 9 (P) = ^Pfc^(p fc ) + inf{Tr[a]} . (26) 

k=0 

Upper bounds on R g (p) can then be obtained by estimating the above minimum value, 
trough specific choices of a. 

A simple possibility for curing the non positivity of ~Pnb is to add to it a matrix pro- 
portional to the modulus of its largest negative eigenvalue; in general, the value of this 
eigenvalue is however difficult to estimate. Another possibility is suggested by the general 

6 This is a direct consequence of the positivity of p, since pb is made of its principal minors. 
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theory of positive matrices (see Theorem 6.1.1 and 6.1.10 in [51]): a sufficient condition for 
a generic hermitian matrix to be positive is that it must be "diagonally dominated", 
i.e. Ma > l-^%l> Then, in a fixed separable basis, by choosing for a the diagonal 

matrix whose entries are given by the sum of the modulus of the elements of the correspond- 
ing rows of — PnB) the matrix a in ( 125)) results positive and makes the combination p + a 
separable. One can then conclude that: 

N 

Rg(p) < ^PkRg(pk) + \ \pNB\\h , (27) 
fc=0 

where, for any matrix M, ||M||^ = ^ . is the so-called £i-norm (see [5Tj)FI 

In presence of just two modes, M = 2, each of which forming a partition, the above 
considerations further simplify. In this case, the Fock basis in (JB]) is given by the set of 
N + 1 vectors {\k; N — k), < k < N}, without the need of further labels; indeed, the N 
bosons can occupy either one of the two modes. Notice that by (j2J) this set of Fock vectors 
constitutes the only basis made of separable pure states [T2] . In this basis, a generic density 
matrix for the system can then be written as: 

N N 

P=^ Pki\k;N -k)(l-N -l\ , ^2p kk = l. (28) 

k,l=0 k=0 

By Proposition 1, once adapted to this simplified case, it follows that a state as in ( 128]) is 
separable if and only if p^i ~ 5m, i-e. the density matrix p is diagonal in the Fock basis. As a 
consequence, an entangled state can never be made separable by mixing it with a separable 
state, so that its robustness of entanglement results always infinite. 

In the case of the generalized robustness, since there are only diagonal and off-diagonal 
terms and no blocks, the above discussed upper bounds (122)) and ( 127)) simplify, becoming: 

i) R g (p)<X{N+l) , ii) R g (p)<\\p ND \\ tl , (29) 

where, as before, pnd is the non diagonal part of p in the Fock basis, and A the modulus of 
its largest positive eigenvalue. These bounds can be explicitly evaluated for specific classes 
of states as shown below. 

Examples: i) Let us first consider pure states of the form: 

1 - 

p=|V)M, W) = -7frr^J2Pk\k-,N-k) , Pk = e^. (30) 

v 7V + 1 k=0 

7 Note that this condition is base-dependent: nonequivalent conditions are obtained by expressing the 
matrix M in different basis. 

8 The same procedure can also be applied to the previously used decomposition of p into its diagonal and 
off-diagonal parts: p = pu + Pnd\ the mixing matrix a would now be composed by — Pnd plus a diagonal 
matrix whose entries are given by the sums of the modulus of the elements of each row of pnd- In this 
case, one easily finds that: R g (p) < | |pjv _d | Ui ; although in general ||/9jvd|Ui ^ ||a/v\b|Uij this constitutes a 
different upper bound for the generalized robustness, independent from that given in (|27l) . 
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The non-diagonal part of the matrix p is given by pnd = ($ — + 1) where $ = 

J2 k ,i e^-w) \k- N - k)(l; N - l\ is the (N + 1) x (N + 1) matrix of phase differences; its 
eigenvalues are zero and N+l, so that the largest negative eigenvalue of — Pnd is in modulus 
N/(N + 1). From the first bound in ( 12"§|) one therefore gets: R g {p) < N. This is also the 
result of the second bound, since the norm HpjvdII^i is also equal to N. 
ii) Nevertheless, the two bounds in (129)) do not in general coincide, as can be seen by 
slightly generalizing the states in (130]) by allowing the coefficients pu to acquire a non unit 
norm, pk = \pk\e llpk , J2k \Pk\ 2 = 1- By choosing the norms and phases of the p^'s to be 
uniformly distributed, one can easily generate states for which the second bound in ( 129|) is 
more stringent. 

Hi) Note, however, that the hierarchy of the bounds in (129 p can be reversed, as it happens 
for instance with the following mixed state: 

N N 

p = y^lfc.N- k)(k; N-k\ -. r > life; iV - k)(l; N — l\ . (31) 

k=0 v 1 k,l=0 

k^l 

Indeed, now one has p^D — (1 — E)/N{N + 1), where E is the matrix with all entries equal 
to one. One easily checks that | |pjvd| Ui — 1> while the modulus of the largest negative 
eigenvalue A of — Pnd is 1/N(N + 1). Therefore, from the first bound in ([29]) one gets 
Rg{p) < 1/-^; which is lower than the second one. □ 



4 On the geometry of iV-boson states 

As shown by the previous results, the properties of the states of a system of iV-bosons result 
rather different and to a certain extent richer then those describing distinguishable particles. 
This is surely a consequence of the indistinguishability of the system elementary constituents, 
but also of the presence of the additional constraint that fixes the total number of particles to 
be N. This additional "rigidity" allows nevertheless a detailed description of the geometric 
structure of the set of iV-boson states. 

Let us first consider the set of entangled states. As mentioned before, for systems of iden- 
tical particles, negativity, as defined in ( Tl0|) . results a much more exhaustive entanglement 
criteria than in systems of distinguishable particles: it seems then appropriate to look for 
states that maximize it. 

Proposition 4. Given any bipartition (m,M — m), the negativity M{p) is maximized by 
pure states that have all the Schmidt coefficients nonvanishing and equal to a normalizing 
constant. 

Proof. First observe that the negativity is a convex function [19] , i. e. it satisfies the inequality 

i i i 
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for any convex combination YliPi Pi °f density matrices p». Since any density matrix can be 
written as a convex combination of projectors over pure states, one can limit the search to 
pure states; in a given (m, M — m)-bipartition, they can be expanded in terms of the Fock 
basis ([HI) as: 

\*P) = ^2\^k), \ipk) = ^ka*>\k,a;N -k,a'} , ^^|^ w | 2 = l. (33) 

fc=0 a=\ a' = l k au' 

As observed before, for each k, the set of vectors {\k,a; N — k, a 1 )} span a subspace Tik C H 
of finite dimension D^- k of the total Hilbert space H, and the component is a vector 
of this space. Recalling Remark 3, one can then write it in Schmidt form, 

\iP k ) =J2^ka\\k,a;N-k,a)) ; (34) 

a=l 

in this decomposition, the orthonormal vectors {\\k, a; N — k, a))} form the Schmidt basis, 
with T>k = min{Dk, Djy^k} (c/. (JED), while the Schmidt coefficients tp^a are non negative 
real numbers, satisfying the normalization condition 'Yl ia {'4 ) ka) 2 — 1- m this new basis, 
one can easily compute the negativity of the density matrix \%p)(%p\ ) 



N V k 
k=0 a=l 



(35) 



Clearly, the negativity increases monotonically with the sum J^fcX/a^fo*' that therefore 
needs to be maximized under the constraint YlkYla^ka) 2 = 1- One easily shows that the 
maximum is obtained when all coefficients ip^a are equal and constant, 

1 N 
i) ka = — = , V = 2J T^k , (36) 

V k=0 

so that, A/'(|V')(' i /'l) = (^ — l)/2- Further, all Schmidt coefficients need to be non vanishing 
in order to get this maximum value for A/"; indeed, for any state with Schimdt number 
&<T>, one has: < (P-l)/2. □ 

In analogy with the standard case, states for which the negativity reaches the maximum 
value (T> — l)/2 can be called "maximally entangled states". 

Remark 5: i) Given a fixed bipartition (m, M — m), let us consider such a maximally 
entangled state p — |^)(^|. By tracing over the degrees of freedom of the modes pertaining 
to the second partition, one obtains the reduced density matrix p^ 1 ' describing the first m 
modes only: 

1 N V k 

p (1) = ^£En*' a » «*' a ih (37) 

k=0 o=l 
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similarly, by tracing over the first partition, one obtains the reduced density matrix p^ 
describing the second N — k modes. In the case of distinguishable particles, either p^ or 
p( 2 ) is proportional to the identity matrix; this is no longer true here: in fact, given a block 
with k fixed, this happens only when < Dn-u-, *.e. when 2\ = D^. 
ii) On the other hand, given the expression ([3"T|) . one easily computes its von Neumann 
entropy, obtaining: S(p^) = — Tr [p^lnp^] = In 22, as in the standard case. 

in) Similarly, the purity of p^ is given by: Tr[(p^^) 2 ] = 1/T>. This result and that of the 
previous remark can equally be taken as alternative, equivalent definitions of the notion of 
maximally entangled states. □ 

Let us now come to the analysis of the space S of separable states as determined by the 
generic bipartition (m, M — m). Among them, the totally mixed state p m i x , proportional to 
the unit matrix, stands out because of its peculiar properties. In fact, recall that in the case 
of distinguishable particles, p mix always lays in the interior of S [HJ [29]; instead, now one 
finds: 

Proposition 5. Given any bipartition (m, M — m) and an associated separable basis made 
of the Fock vectors \k,a; N — k, a'), the totally mixed state, 

1 - 

Pmix = -p^ Yl \k,a;N-k,a')(k,a;N-k,a'\ , (38) 
lays on the border of the space S of separable states. 

Proof. Let us take a state p en t which can not be written in block diagonal form as in ffTTj) . 
f lT2|) : by Proposition 1, it is entangled. Then, consider the combination 



Pe = Y^(/°mix + ePent) , 6 > . (39) 

For any e, this combination will never be separable, since it is not block diagonal. In other 
terms, in the vicinity of p m ; x , one always find entangled states, so that it must lay on the 
border of S. □ 

Remark 6: i) Note that similar considerations apply to all separable states: there always 
exist small perturbations of separable, necessarily block diagonal, states that make them not 
block diagonal, hence entangled. Instead, in the case of distinguishable particles, almost all 
separable states remain separable under sufficiently small arbitrary perturbations [29, 152]. 
ii) Further, using analogous steps, one can show that a Werner-like state, pw = P |V')(V ; I + (1 — 
p) Pmix, < p < 1, with a maximally entangled state, is entangled for any nonvanishing 
value of the parameter p, while for distinguishable c?-level particles, this happens only when 
l/(d+ 1) <p < 1 [53]. □ 

The result of Proposition 5 is most strikingly illustrated by considering a system of iV 
bosons that can occupy two modes (M = 2), each of which forming a partition. In the 
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previous section, we have seen that a generic state p for this system can be written in terms 
the Fock basis as in f )28|) ; further, it results separable if and only the density matrix p is 
diagonal in this basis. Instead, p € given in fl39|) will develop non-diagonal entries as long as 
e starts to be non vanishing. 

As mentioned before, in this case the set of N + 1 vectors {\k; N — k)} constitutes the only 
basis made of separable pure states [12]. Therefore, the decomposition of a generic separable 
state p in terms of projections on separable states turns out to be unique. As a consequence, 
the set S of separable states is a sub-variety of the convex space M. of all states that turns 
out to be a simplex, whose vertices are given precisely by these projections. 

The space of states of N bosons confined in two modes is then much more geometrically 
structured than in the case of systems made of distinguishable particles. Indeed, given a 
complete set of observables {Oi}, one can decompose any density matrix p as: 



P 



this decomposition is over pure states, whenever the operators Oi are chosen to be projectors. 
Therefore, in general, there are infinite ways of expressing a density matrix as a convex 
combination of projectors, even when these projectors are made of separable states. As seen, 
this conclusion no longer holds for systems of identical particles. 
The totally mixed state p mix has also another interesting property^ 

Proposition 6. The totally mixed state is the only state that remains separable for any 
choice of bipartition. 

Proof. First, let us again consider the simplest case M = 2. Any Bogolubov transformation 
maps the set of creation and annihilation operators a\ and cij, % = 1,2 into new ones b\, 6j, 
i = 1, 2; a simple example is given by 

a 1 + a 2 a l -a 2 , , 

= — • h = — • (41) 

and their hermitian conjugates. The operators b\, b; L define a new bipartition (B±, B 2 ) of the 
full algebra B of bounded operators, distinct from original one (Ai,A 2 ) generated by a\, a t . 
States (and operators as well) that are local in one bipartition might turn out to be non-local 
in the other. For instance, the Fock states {\k, N — k)} result entangled with respect to the 
new bipartition defined by the transformation f )4T|) ; indeed, one finds: 

!*■"-*> = ^ J. W EE (!) I"' (4f — 10) , (42) 



2«ft ^k\(N - *)! ^ ^ V 



9 For systems of distinguishable particles, the problem of finding so-called "absolutely separable states" , 
i.e. states that are separable in any choice of tensor product structure, has been discussed in [511122]. 
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so that \k,N — k) is a combination of (Bi, i32)-separable states. A similar conclusion applies 
to the mixed states (125]): in general, any separable state p sep = pk \ k; N — k) (k; N — k\ is 
mapped by a Bogolubov transformation into a non-diagonal density matrix, and therefore 
into an entangled state. In fact, one can always find a unitary transformation U that maps 
any diagonal matrix p into a non-diagonal one UpU^ . In particular, when p is not degenerate, 
the transformed matrix UpU^ results diagonal only if the operator U is itself a diagonal 
matrix; in this case however the corresponding Bogolubov transformation results trivial and 
does not define a new partition. The only density matrix that remains invariant under all 
unitary transformations is the one proportional to the unit matrix, i.e. the totally mixed 
one. This conclusion can easily be extended to the multimode case: given any separable state 
in a given (m, M - m) bipartition, p sep = Y,k=o J2a,a' Pkaa>\k, a;N - k, a')(k, a; N - k,a'\, 
one can always construct a Bogolubov transformation that maps it into an entangled one: it 
is sufficient to apply the above considerations to any couple of modes belonging to separate 
partitions. Therefore, also in this more general setting, only the state p mix in ( 155]) is left 
invariant by all Bogolubov transformations. □ 

Thanks to these results, the global geometry of the space of iV-boson states starts to 
emerge more clearly. Again the two-mode case is easier to describe. We have seen that by 
fixing a bipartition one selects the sets S of separable states; these form a sub-variety of the 
convex space Ai of all states, forming a (N + l)-dimensional simplex, with the projectors 
over Fock states as generators. Changing the bipartition through a Bogolubov transformation 
produces a new simplex, having only one point in common with the starting one, the state 
p mix . The geometry of the space of two-mode iV-boson states has therefore a sort of star-like 
topology, with the various simplexes sharing just one point, the totally mixed state. 

The case of M modes is more complex. For a fixed bipartition, the space of separable 
states is a sub-space of the convex space of all states which is not any more strictly a simplex: 
the decomposition of generic separable state p is no longer unique, since the Fock states in 
([6]) are no longer the only separable pure states: for each k, reshufffings over the indices 
a, a' are still allowed. Nevertheless, also in this case the global state space presents a sort 
of star-like topology: only one point is shared by all separable bipartition sub-spaces, the 
totally mixed state p mix . 

5 Outlook 

In many-body systems composed by a fixed number of identical particles, the associated 
Hilbert space no longer exhibits the familiar particle tensor product structure. The usual 
notions of separability and entanglement based on such structure is therefore inapplicable: 
a generalized definition of separability is needed and can be given in terms of commuting 
algebras of observables instead of focusing on the particle tensor decomposition of states. 
The selection of these algebras is largely arbitrary, making it apparent that in systems of 
identical particles entanglement is not an absolute notion: it is strictly bounded to specific 
choices of commuting sets of system observables. 
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Using these generalized definitions, we have studied bipartite entanglement in systems 
composed by N identical bosons that can occupy M different modes. More specifically, we 
have analyzed to what extent entangled states result robust against mixing with another 
state, either separable or entangled. We have found that in general, the entanglement con- 
tained in bosonic states is much more robust than the one found in systems of distinguishable 
particles. This result has been obtained by analyzing the so-called robustness and general- 
ized robustness of entanglement, of which explicit expressions and upper bounds have been 
respectively given. A quite general characterization of the geometry of the space of bosonic 
states has also been obtained: this space exhibits a star-like structure composed by intersect- 
ing subspaces each of one determined by a given bipartition through the subset of separable 
states. All these separable subspaces share one and only one point, the totally mixed state, 
hence the star-shape topology. 

As a final remark, notice that all above results can be generalized to the case of systems 
where the total number of particles is not fixed, but commutes with all physical observables 
(i.e. we are in presence of a superselection rule [55]). In such a situation, a general density 
matrix p can be written as an incoherent mixture of states pn with fixed number N of 
particles, having the general form fl9]); explicitly: 

P = ^2 ^ nPn ' >0 , ^ Aat = 1 . (43) 

N N 

The state p is thus a convex combination of matrices pn having support on orthogonal 
spaces. Arguments similar to those used in proving Proposition 2 (and Proposition 3) allow 
us to conclude that both the robustness and the generalized robustness of entanglement 
of the state p in fl4"3"j) is the weighted average of the robustness of the components p^, i-e. 
R(p) = J2n R(pn) f° r both cases. The problem of computing the robustness of incoherent 
particle number mixtures (1431) is then reduced to that of determining the robustness of the 
corresponding components at fixed particle number, for which the considerations and results 
discussed in the previous sections apply. 
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